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Preface 

This  analysis  of  the  M61  movable  gun  is  intended  to  be  an  overall 
evaluation  of  the  control  system.  Originally,  the  concept  of  the  thesis 
was  to  design  a  digital  controller  for  the  system.  But  the  gun  servo 
subsystem  dynamics  are  good  enough  and  the  sampling  rate  so  low  that  a 
conventional  digital  controller  is  neither  required  nor  feasible.  Thus 
the  motivation  became  one  of  answering  questions  which  were  raised  in 
the  review  of  readily  available  literature. 

The  gun  system  analysis  is  based  on  locally  available  Delco  Cor¬ 
poration  and  McDonnell  Aircraft  Company  reports.  Mo  attempt  was  made 
to  contact  either  contractor  regarding  specific  questions,  so  portions 
of  the  given  information  may  not  be  in  agreement  with  current  data  or 
design. 

1  extend  my  appreciation  to  my  thesis  advisor,  Capt  James  Silver- 
thorn© ,  and  Professor  C.  H.  Houpis  for  their  assistance  throughout  this 
study  and  to  Professors  D.  W.  Breuer  and  Robert  Calico  for  their  interest 
in  reviewing  my  thesis.  I  am  also  indebted  to  Lt  Col  tathony  Leathern 
who  sponsored  this  thesis  and  to  Joe  Rogers  and  Major  K.  E.  Hudson  who 
provided  additional  background  material. 

Finally,  and  most  dearly,  I  would  like  to  thank  my  wife,  Pat,  for 
the  time  she  spent  with  and  without  me  over  the  course  of  my  studies. 

Her  encouragement  and  understanding  were  invaluable. 


Donald  E.  Jones 


ii 


Contents 

Page 

Preface  .  .  . . ii 

List  of  Figures  . . . . . .  v 

List  of  Tables  .  .................  vli 

Symbols . viii 

Abstract . x 

I  Introduction  ......  .  ......  .  1 

Historical  Background  ...  .  ..........  1 

System  Definition . ••••• .  2 

Purpose  . . 2 

Assumptions . .  4 

Performance  Criteria  . . 5 

Approach  .  ..........  .  7 

II  Technical  Description  of  the  Gun  System .  8 

Gun  Servo  Subsystem  . . 10 

Digital  Controller  . . 16 

Coordinate  Transformation  .  ...........  17 

III  Analysis  of  the  Gun  Servo  Subsystem  . . 21 

Differential  Pressure  Compensator  .  21 

Feedforward  Compensator  .....  .  .......  30 

Effect  of  Rate  Conmands  and  Rate  Feedforward . 34 

IV  Analysis  of  Digital  Portion  of  Gun  System . .36 

System  Gain  .  .................  36 

Digital  Rate  Feedforward  42 

Sensor  Filters . 42 

Computation  Time  ....••••• . 46 

Word  Length  Considerations  ••.•••.••••••••  31 

V  Analysis  of  Muzzle  Response  ••••••••••.••...  32 

Effect  of  Structural  Modes  on  Muzzle  Response  .  32 

Controller  for  Structural  Modes  ...••••.•••..  36 

VI  Conclusion . 66 


Suonary  •  .  .  • 
Recommendations 


66 

67 


9  *«g« 

Bibliography  ..............  .  .  .  68 

Append  ixA:  Derivation  of  Gun  Dynamics  Model  . .  70 

Appendix  Bt  Determination  of  Gun  Output  Coefficients . .  73 

Appendix  C:  Gun  Servo  Subsystem  State  Variable  Representation  •  •  81 
Appendix  Di  M6i  Simulation  Program  ••..•••••••••••  88 

Appendix  Es  Z  and  S  Plane  Relationships 


108 


Cigna 


Cam 

3 


1  M61  Movable  Gun  Servo  Subsystem 


2  Basle  Gun  System  Model 


3  Gn  System  Model 


4  Gun  Servo  Subsystem  Model 
3  Sensor  Frequency  Response 


6  Gun  Geometry 


7  Gun  System  Interface 


8  Continuous  Complementary  Filter 


9  Baseline  Gun  System  Response 


10  Differential  Pressure  Compensation  Loop 


11  Differential  Pressure  Compensator  Root  Locus 


12  Differential  Pressure  Loop  Time  Responses 


13  Gun  Servo  Subsystem  Reduction 


14  Gun  Servo  Subsystem  Root  Locus 


IS  Gun  Servo  Subsystem  Response  Using  Position  Inputs 


16  Discrete  Gun  System  Model 


17  Discrete  Gun  System  Root  Locus 


18  Response  of  Gun  System  with  No  Digital  Rate  Feedforward  ,  .  . 


19  Effect  of  10  Hz  Sensor  Filter  an  System  Response 


20  Gun  System  Model  Including  Computational  Delay 


21  Effect  of  Computation  Time  on  System  Response 


22  Gun  Structural  Mode  Slopes 

23  Muzzle  Frequency  Response  . 


24  Muzzle  Time  Rezponse 


25  Structural  Mode  Compensator 


59 


Figure  Page 

26  Structural  Mode  Compensator  Root  Locus .  60 

27  Structural  Mode  Compensator  Discrete  Root  Locus  ••••••  61 

28  Structural  Mode  Compensator  Time  Response .  62 

29  Modeled  Muzzle  Frequency  Response  .  ••••••••  74 

30  Modeled  Sensor  Frequency  Response  . . 79 

31  State  Variable  Block  Diagram  81 

32  Z  Domain  Unit  Circle  . . .  •  •  •  •  108 


I 


I 


5 


List  of  Tables 

able  Pag 

Gtm  Servo  Subsystem  Specifications . .  .  *  •  . 

II  Differential  Pressure  Compensator  Gin  Response  .•••.••  25 

III  Differential  Pressure  Compensator  System  Response  .  .  •  •  •  25 

IV  Gun  Servo  Subsystem  Ramp  Errors  33 

V  Effect  of  Kql  on  Gun  System  Performance . 33 

VI  Effect  of  Rate  Feedforward  on  Gun  Response . 34 

VII  Comparison  of  Analytical  and  Simulation  Response  ......  41 

VIII  Poles  of  Closed  Loop  Gun  System  as  a  Function  of  Gain  ...  41 

IX  Effect  of  Sensor  Filter  on  System  Response  .........  44 

X  Effect  of  Digital  Filter  Gain  on  System  Response . 46 

XI  Comparison  of  Simulation  and  Analytical  Computation 

Time  Effects  .........  .  .........  49 

XII  System  Response  for  Different  Computation  Times . 51 

XIII  Structural  Mode  Compensator  Response . 64 


Tii 


Svnbels 


Gun  dynamics  output  coefficient  for  1th  mode 
Computation  time  delay  (Tc) 

Actuator  piston  force 
Sampling  frequency  (Herts) 

Feedforward  compensation 
Servovalve  dynamics 
Gun  dynamics 

2  mode  gun  dynamics  model 
4  mode  gtn  dynamics  model 
Muzzle  dynamics  model 
Differential  pressure  compensator 
Sensor  filter  (analog) 

Gun  system  transfer  function 

Gun  servo  subsystem  transfer  function 


Gun  dynamics  with  differential  pressure  compensator 

Zero  order  hold  dynamics 

Differential  pressure  compensator 

Feedback  structural  compensator 

Sample  period 

System  gain 

Digital  sensor  filter  gain 
Gun  servo  subsystem  gain 
X  actuator  length 
Y  actuator  length 


Till 


Centered  x  actuator  length 


•  So 

So 

*p 

■rad 

MSE 

*P 


*c 

b 

*s 


x 


*s 

*M 

*C 

*Cd 

T. 

*M 

s 

Z 

e 

®c 

» 


» 


Centered  y  actuator  length 
Peak  overshoot 
Milliradlans 
Mean  Squared  Error 

Input  to  compensated  servovalve-actuator»gun  subsystem 

Laplace  domain  variable 

Sample  time 

Computation  time 

Peak  overshoot  time 

Rise  time 

Settling  time 

State  variable;  X  actuator  position 
Estimated  position*  Estimated  state 
Sensor  output  (X  channel) 

Muzzle  output  (X  channel) 

GSS  rate  command 
Digital  rate  command 
Sensor  output  (Y  channel) 

Muzzle  output  (Y  channel) 

Discrete  domain  variable 
Z  transform  operator 
Elevation 
Elevation  command 
Natural  frequency 
Azimuth 

Azimuth  command 


♦ 

♦c 


lx 


{  Damping  ratio 

A  1  -  C/T 


AFIT/GA/AA/78D-5 


Attract 


'  A* 

/  ” 


The  effects  of  changing  control  parameters  of  the  movable  M61  gun 
system  proposed  for  the  F-15  aircraft  are  examined  using  time  response 
and  root  locus  methods.  In  the  course  of  the  analysis,  a  Fortran  IV 
simulation  program,  state  space  model,  and  gun  servo  subsystem  Z  trans¬ 
form  are  developed. 

The  gun  servo  subsystem  design  has  little  effect  on  system  response. 
The  system  settled  in  under  0.2  sec  and  had  less  than  10%  overshoot  for 
any  open  loop  gain  from  0  to  200  sec"1  and  with  or  without  differential 
pressure  compensation. 

The  overall  system  is  stable  for  a  system  gain  of  0  to  39  and 
exhibits  nearly  deadbeat  responses  for  a  gain  of  20.  Digital  rate  feed¬ 
forward  is  required  to  keep  ramp  following  error  below  1  mrad  for  a 
5°/sec  ramp.  Digital  filtering  improves  response  and  analog  low-pass 
sensor  filters  with  a  cutoff  of  30  Hz  eliminate  aliasing  while  moderately 
reducing  system  performance.  Computation  delays  of  less  than  0.005  sec 
were  found  to  have  negligible  effect  on  the  system  response. 

The  muzzle  response  is  examined  and  a  compensator,  which  neglects 
barrel  cluster  rotation, is  designed  to  reduce  the  50%  overshoot  and  over 
2  sec  settling  time  for  a  step  input.  This,  however,  degraded  tracking 
of  more  realistic  (lower  frequency  content)  inputs  indicating  that  a 
better  compensator  should  be  designed  or  that  muzzle  response  at  target 
acquisition  should  be  allowed  to  settle  before  firing. 

Overall,  the  movable  M61  was  found  to  be  an  extremely  fast  gun 
system,  insensitive  to  most  control  parameters. 
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ANALYSIS  OF  A  CONTROLLER  FOR  THE 


M61  MOVABLE  GUN 


I  Introduction 


Historical  Background 

Since  World  War  I*  fighter  aircraft  have  been  armed  with  rapid 
firing  guns  to  perform  the  close  in  air  superiority  role.  Although  it 
has  been  suggested  that  air  superiority  aircraft  can  function  without  a 
gun,  experience  has  shown  that  the  gun  has  a  place  on  the  highest  tech¬ 
nology  aircraft  (Ref  2:1).  Gunsights  have  been  Improved  dramatically 
since  the  use  of  cross  hairs  and  the  gun  Itself  has  been  improved  over 
the  years.  However,  the  same  method  of  aiming  the  gun  is  still  being 
used,  i.e.,  point  the  aircraft. 

Although  this  has  proven  effective  in  the  past,  increasing  speed 
and  maneuverability  have  put  an  extremely  heavy  burden  on  the  pilot  to 
track  the  target  with  his  aircraft.  If  a  means  were  available  to  relieve 
the  pilot  of  a  portion  of  this  task  and  avoid  the  dynamic  constraints  of 
the  aircraft,  aircraft  gunnery  could  be  much  more  effective. 

The  movable  gun  concept  was  investigated  by  the  RAND  Corporation 
In  1968  (Ref  16).  This  study  indicated  that  the  movable  gun  greatly 
improved  firing  opportunities.  In  addition,  the  greatest  performance 
increase  occurred  within  the  first  few  degrees  of  movement. 

The  movable  gun  concept  was  further  explored  in  the  EXPO  series  of 
air  to  air  fire  control  studies  performed  by  the  McDonnell  Aircraft 
Company  (MCAIR).  EXPO  V  not  only  confirmed  the  RAND  findings,  but 
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generated  a  preliminary  hardware  design  for  a  movable  M61  cannon  for  tho 
F-15  aircraft. 

System  Definition 

The  movable  M61  was  designed  by  Delco  Electronics  Division  to 
MCAIR  specifications.  The  Gun-Servo  Subsystem  (GSS),  consisting  of  the 
gun,  hydraulic  actuators  and  associated  hardware  and  electronics,  is 
shown  in  Figure  1.  The  gun  servo  subsystem  Includes  a  failure  monitor 
which  centers  the  gun  if  an  error  is  detected  in  its  response.  This 
portion  of  the  gun  servo  subsystem  will  not  be  included  in  any  of  the 
analysis. 

The  gun  system  is  composed  of  the  gun  servo  subsystem,  portions  of 
the  F-15  mission  computer  (MC)  and  the  connecting  data  bus. 

Purpose 

The  purpose  of  this  study  is  to  examine  the  effects  of  varying  con¬ 
trol  parameters  of  the  gun  servo  subsystem  as  designed  by  Delco,  examine 
digital  components  of  the  system,  and  examine  time  response  of  the  muzzle. 
Delco  dealt  primarily  with  frequency  domain  responses  in  their  reports; 
this  study  will  focus  on  time  domain  effects. 

At  this  time  there  have  been  no  final  descriptions  of  a  digital 
controller  algorithm  appearing  in  Delco  or  MCAIR  literature.  This  study 
will  examine  some  of  the  aspects  of  the  digital  controller  Including 
open  loop  gain,  digital  rate  feedforward  and  digital  filtering  of  sen¬ 
sor  outputs. 

Finally,  the  effects  of  the  gin  structural  modes  an  the  muzzle 
response  will  be  examined.  The  muzzle  frequency  response  is  given  by 
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Delco  reports  but  Its  time  domain  consequences  are  not  discussed*  A 
preliminary  design  for  a  compensator  to  reduce  structural  resonance 
will  be  discussed. 


Assumptions 

Ihls  thesis  will  sake  several  simplifying  assumptions  regarding  the 
gun  dynamics  and  command  Inputs.  These  assumptions  and  their  implica¬ 
tions  follow. 

Gun  Dynamics.  The  gun  will  be  treated  as  a  single  body,  neglecting 
barrel  cluster  rotation  and  projectile  motion.  The  six  barrel  gun  has 
a  firing  rate  of  6000  rounds  per  minute,  so  gun  firing  occurs  at  100  Hz 
and  the  barrels  rotate  at  16.7  revolutions  per  second.  Gun  firing  will 
be  an  impulsive  input  along  the  gur.  axis,  and  therefore  should  have  a 
negligible  effect  on  the  system  response. 

Since  the  barrel  cluster  rotation  is  a  relatively  low  frequency 
effect  (the  first  structural  mode  is  at  11.6  Hz),  it  may  cause  some 
coupling  of  the  gun  dynamics  in  the  azimuth  and  elevation  axes.  This 
effect  will  have  its  greatest  impact  on  the  structural  mode  analysis. 

By  neglecting  nonlinear  friction  and  treating  structural  damping 
in  the  gun  as  viscous  damping,  the  system  can  be  examined  using  linear 
models.  The  use  of  linear  models  greatly  simplifies  the  analysis. 

Command  inputs.  It  is  assumed  that  rate  and  position  commands  for 
the  gun  azimuth  and  elevation  are  available  within  the  mission  computer. 
It  is  not  the  intent  of  this  study  to  become  involved  with  processes  of 
target  and  projectile  prediction  or  aircraft  dynamics. 


The  Delco  reports  pieced  emphasis  on  frequency  domain  criteria  for 
the  gun  servo  subsystem.  This  thesis  will  deal  almost  exclusively  with 
time  domain  analysis.  Most  of  the  specifications  supplied  by  MCAIR  are 
related  to  gun  servo  subsystem  performance  rather  than  to  the  gun  system 
which  contains  the  digital  control  loop.  Time  domain  specif icatlons 
must  then  be  developed  for  the  entire  gun  system  based  on  the  gun  servo 
subsystem  requirements  and  same  engineering  judgment. 

The  available  performance  specifications  are  listed  in  Table  I. 

When  a  specification  appears  in  more  than  one  reference*  the  primary 
source  is  cited*  When  in  disagreement*  the  most  recent  source  is  U3ed. 


Table  I 

Gun  Servo  Subsystem  Specifications 

Specification _ Reference 

Steady  State  Oscillation  <  0,1%  14 

GSS  Overshoot  <  20%  3 

Bandwidth  of  GSS  100  He  5 

Static  Accuracy  .3  mrad  14 

Rate  Following  Error  (GSS)  <  1  mrad  for  5°/sec  3 

Angular  Excursion  ±  3°  14 

Angular  Velocity  >  4S°/sec  14 

Angular  Acceleration  >  200  Rad/sec^  3 

80%  of  100  Round  burst  in  8  mrad  Dia  Circle  13 


✓ 


The  steady  state  oscillation  and  static  accuracy  arc  primarily 
related  to  hardware  tolerance  and  saisor  errors,  so  will  not  be  used 
except  In  a  short  discussion  of  digital  word  length. 

The  GSS  bandwidth,  overshoot,  and  rate  following  error  specifica¬ 
tions  will  be  addressed  when  discussing  the  Gun  Servo  Subsystem.  The 
bandwidth  of  100  Hs  cannot  be  applied  to  the  discrete  system  since  the 
sampling  rate  Is  20  Hs. 

The  GSS  overshoot  criterion  seems  to  be  high  for  the  system  over¬ 
shoot  based  on  the  dispersion  specification,  so  will  not  be  used  In 
system  analysis.  The  rata  following  specification  will  be  used,  however. 

The  angular  excursion,  velocity,  acceleration  are  functions  of  the 
hydraulics  and  actuators  so  are  not  used  as  control  criteria,  but  as 
limits  In  the  gun  simulation. 

From  the  dispersion  specification,  the  maximum  error  which  will 
place  the  target  In  the  area  of  a  probable  hit  Is  4  mrad.  Based  on  this, 
it  would  be  desirable  to  keep  the  overshoot  less  than  4  mrad  to  maximise 
probability  of  a  hit.  For  the  maximum  excursion  of  52  mrad  this  is  about 
10X.  Since  the  gun  has  a  dispersion  of  4  mrad,  a  5*  settling  criterion 
(2.5  mrad  at  maximum  excursion)  is  used  rather  than  2%.  It  would  be 
desirable  to  have  the  gun  settle  as  fast  as  possible  but  0.2  seconds 
(20  rounds)  sesms  reasonable. 

The  design  criteria  to  be  used  for  the  system  are  than* 

T,  (5%)  0.2  sec 

Mp  <10* 

Kate  following  error  £  1.5* 
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Approach 

Tha  GSS  and  gut  system  will  be  analysed  using  both  analytical 
methods  and  a  computer  simulation.  For  the  analytical  portion  of  the 
analysis,  the  interactive  computer  aided  design  program  TOTAL  (Ref  10) 
vas  used.  The  simulation  program  (described  in  Appendix  0)  is  a 
FORTRAN  IV  program  executed  an  the  CDC  6600/CYBER  digital  computer. 

The  computer  simulation  includes  some  of  the  nonlinearities  of  the  sys¬ 
tem  and  uses  more  complete  gun  dynamics. 
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II  Technical  Description  of  the  Gun  System 

The  basic  gin  system  can  be  modeled  as  shown  in  Figure  2.  This  non* 
linear  servo  system  is  driven  by  azimuth  and  elevation  cousands  internal 
to  the  F-15  mission  computer.  These  commands  are  generated  using  outputs 
of  the  lead  computation  routines  within  the  mission  computer  along  with 
target  tracking  information  from  the  AFG-63  Radar.  As  a  result  of  the 
limited  data  rates  of  these  inputs  and  the  computational  burden  of  the 
mission  computer*  position  coomands  are  available  at  a  20  Hz  rate. 

A  more  detailed  gun  system  model  is  shown  in  Figure  3.  In  the 
following  sections*  this  model  will  be  broken  down  and  each  of  its  com¬ 
ponents  described. 


cmd  '< 

MISSION 


COMPUTER 


t'cmd 


X  ACTUATOR 


i  Y  ACTUATOR 


GUN 

GEOMETRY 


0 


digital 

CONTROL 


GUN  SERVO  SUBSYSTEM 


Figure  2*  Basic  Gin  System  Model. 
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Figure  3.  Gun  System  Model 


Gun  Servo  Subsystem 


The  gun  servo  subsystem  contains  two  nearly  orthogonal  actuator 
channels*  Each  channel  consists  of  the  actuator,  servovalve,  sensors 
and  the  associated  compensation  networks.  A  model  for  each  of  the 
identical  channels,  including  gun  dynamics,  is  shown  in  Figure  4  (Ref  3). 


Figure  4.  Gun  Servo  Subsystem  Model. 


Feedforward  Compensation.  The  feedforward  compensation,  Gj^(s),  Is 
simply  a  gain  in  the  current  design* 

G|(s)  -  K(,L  -  100  sec*1  (1) 

Other  compensators  and  gains  will  be  discussed  in  chapter  III. 

Servovalve  Dynamics.  The  servovalve  dynamics  relate  electrical 
inputs  to  a  force  exerted  on  the  actuator  piston.  A  model  of  these 
dynamics  is  given  by  the  manufacturer  as 


C2<s) 


(2) 


Ay  (l)y^ 

8(S2  ♦  2  ^  (dyS  ♦  (Ily2) 

where 

Ay  -  l 

ly  “  .70 

(dy  -  200  Hz  ■  1256  rad/sec 

CtTl  pypapi<<.s-  The  gun  dynamics  represent  the  relationship  between 
the  actuator  forces  and  the  gun  displacement.  Because  the  gun  is  not  a 
rigid  body,  a  finite  element  analysis  was  performed  by  Delco  to  obtain 
the  structural  response  of  the  gun  and  actuator  body.  The  finite  ele¬ 
ment  program  provided  the  parameters  for  use  in  an  elaborate  gun  model 
described  in  Ref  13.  The  program  also  produced  the  frequency  response 
of  the  gun  sensor  to  an  actuator  input  shown  in  Figure  5.  This  fre¬ 
quency  response  is  used  to  generate  a  model  for  the  gun  dynamics. 

Two  sensors  are  being  considered  for  the  gw  system.  The  first  is 
a  linear  variable  differential  transformer  (LVDT)  which  senses  actuator 
length.  The  second  is  an  angular  resolver  which  measures  gun  angle  at 
the  pivot  position.  The  gw  has  nearly  Identical  frequency  respwse  at 
both  of  these  sensor  locations  so  the  IV DT  respwse  can  be  used  for  both 
sensors. 

Two  sets  of  gw  dynamics  based  w  this  frequency  respwse  were  used 
in  this  gun  system  analysis.  The  first  is  a  model  used  by  Delco  in 
their  gun  servo  subsystem  design  and  is  used  in  the  analytical  analysis. 
The  second  is  a  more  detailed  model,  using  the  first  four  structural 
modes,  and  is  used  in  the  computer  simulatiw. 
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where 

(j)j  m  ll.S  He  -  72.26  rad/sec 
(t>2  ■  95  He  -  596.9  rad/sec 

(i)_  -  12  He  »  75.4  rad/sec 

U 

t  -  .15 

A  more  detailed  gun  model  is  obtained  by  treating  the  gun  as  a 
structurally  damped  beam.  There  are  many  ways  of  modeling  structural 
damping.  While  viscous  damping  is  one  of  the  least  accurate  for  a  steel 
structure,  it  does  allow  the  use  of  a  linear  model. 

Equation  (4)  (derived  in  Appendix  A)  models  the  gun  dynamics  as  a 
sum  of  Mn"  second  order  modes. 


G3  (s) 


*SS±  .  5  _ Al 

f(s)  i»i  s2  *  2l  ♦  (t>i2 


(4) 


where  (dj  is  the  frequency  of  the  ith  mode*  A^  is  the  output  coefficient 
of  the  1th  mode  and  /  the  gun  damping  ratio. 

The  first  four  natural  modes  (n<»4)  of  the  gun  were  used  in  the 
computer  model.  The  frequencies  of  these  modes  are  given  as  (Ref  5): 

W*  •  11.6  Hz  -  72. 88  rad/sec 

(i)2  ■  34.9  Hz  ■  219.7  rad/sec 

(1)3  -  68.7  Hz  »  431.6  rad/sec 

(1)4  <■  83.3  Hz  m  524.6  rad/sec 

The  output  coefficients  (Ai  -  A^)  and  damping  ratio  (£  )  are  found 
in  Appendix  B  using  the  gun  frequency  response.  These  coefficients 
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are  I 


Aj  —  356*2 
A2  »*  0 

A3  -  48170 
A^  -  185600 

{  -  0.05 

Differential  Pressure  Compensator.  The  differential  pressure  com¬ 
pensator  was  developed  by  Delco  to  provide  damping  of  the  68  and  83 
hertz  gun  structural  modes.  This  compensator  uses  as  its  input  the 
pressure  differential  across  the  actuator  piston.  This  compensator  is 
given  by 

Hpfe)  .  ^g4(s)  (5) 

A 

where  Khyd  is  the  hydraulic  spring  constant  and  A  is  the  piston  area. 
Kfcyd  .  100,000  lb/  in 

A  .9  in^ 


where 


G4(s) 


kAp  <Tj) 

(TjS  ♦  1)(t2s  ♦  1) 


m  3.15  x  10“^  ln/sec/ps i 
l/T\  -  500 
l/T2  -  2500 


(6) 


Gun  Geometry.  The  actuators  are  mounted  on  the  rear  of  the  M61  gun 
as  shown  in  Figure  6.  The  transformation  from  actuator  lengths  to  gun 
angles  is  given  by  (Ref  5) 
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a 

d 


sin  Tjf  -  ~  [2I*2  ly2*  2(2b)\Z+  2(2b)2I*2-Z*4-Iy4-(2b)4]%. 

where  0  is  the  elevation  angle,  tj;  the  azimuth  angle,  1^  and  ly  the 
x  and  y  actuator  lengths,  respectively,  and 
a  ■  b  »  10.25  in 
d  *•  24.1  in 

While  these  transformations  are  not  exact,  they  have  a  maximum  error  of 
0.002  mrad  over  the  3°  excursion. 


Figure  6.  Gun  Geometry 


Interfaces  to  Mission  Control  Computer.  The  gun  servo  subsystem 
also  includes  interface  elements  to  the  mission  computer.  Those  of 
interest  in  the  control  analysis  are  shown  in  Figure  7.  The  interface 
contains  10  bit  analog  to  digital  and  digital  to  analog  converters,  a 


RATE 

COMMAND 


ZOH 

-  ^ 

/ 

SMOOTHING 

FILTER 


6gss(s) 


I 


SENSOR 

t 

GUN 

FILTER 

GEOMETRY 

INTEGRATOR  OUTPUT 


MEASURED 
GUN  AN'GlE 


ZOH  t  Zero  Order  Hold 
G  i  Gun  Servo  Subsystem 

gSo 


Figure  7.  Gun  System  Interface. 

zero  order  hold,  and  two  low  pass  filters.  One  of  these  filters  is 
used  to  smooth  zero  order  hold  output  and  the  other  is  a  filter  for  the 
sensor.  The  smoothing  filter  was  not  included  in  the  simulation  due  to 
an  oversight.  However,  same  of  the  filter  effect  can  be  determined 
based  on  analysis  of  the  sensor  filter  since  the  two  are  essentially  in 
series.  The  integrator  output  is  fed  back  to  the  mission  computer  for 
a  digital  filtering  algorithm. 


The  digital  portion  of  the  gun  system  was  not  well  defined.  Figure 
3  contains  a  composite  of  the  controllers  described  in  References  5,  12, 
and  13.  The  location  of  the  coordinate  transformations  is  arbitrary; 
however,  the  specific  transform  used  is  dependent  upon  their  location. 


16 


f. 

k  ! 

r  ' 

fc  I 

>  I 


I 


The  system  gain  is  given  by  Kj.  The  digital  rate  feedforward  is 
included  in  the  system  per  Reference  12,  This  rate  information  could  be 
generated  inside  the  controller  by  differentiating  the  position  input 
signal,  but  should  be  available  from  the  tracker/predictor  algorithm. 

The  actuator  conutands  are  limited  to  ±  18  inches/sec  corresponding  to 
a  45°/sec  angular  rate. 

Coordinate  Transformation 

There  are  essentially  two  coordinate  transformations  employed  in 
this  version  of  the  digital  controller.  The  first  is  for  the  position 
commands  and  measurements.  The  actuator  positions  are  given  by 

*  -  **  “  **o 

(8) 

Y  "  b  m  V 

where  I*0  and  Ly0  aw  the  neutral  position  actuator  lengths  and  I*  and 
1^  are  the  current  lengths. 

*xo  ■  I^o  "  J  sf  +  b2  <9) 


The  actuator  lengths  are  given  by  (Ref  5) 

Ijc  ■  C  (ty  8  +  ♦  (®  d  •  b)^ 

1^  •  C  CJj  d  ♦  a)^  ♦  (0  d  ♦  b)^ 


(10) 


where  a,  b,  and  d  are  as  defined  above.  Note  that  this  transformation 
uses  the  small  angle  approximations  for  the  sine  so  is  less  accurate 
than  the  gun  geometry  transformation  of  £q  (7). 
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/ 


* 


X(s) 


— i-  Xg(s)  ♦  -a-  xx(8) 
T*+l  T8+1 


From  this*  it  can  ba  seen  that  the  filter  generates  an  estimate  of  the 
gun  position  by  summing  the  low  frequency  components  of  the  sensor 
measurements  and  the  high  frequency  components  of  the  integrator  output. 
A  discrete  version  of  Eq  (12)  must  be  developed  for  implementation  in 
the  digital  computer. 

Since  the  inputs  to  the  digital  filter  are  discrete  Impulses,  a 
sero  order  hold  must  be  placed  an  the  input  to  D(s).  Then 

D*(s)  -  gzqh(s)  D(s) 


X*(s)  -  0<s)Xa*<»)  -  D’(s)Xi*(s)]  ♦  Xi(s)  (13) 

where  the  star  indicates  a  sampled  input.  Since  the  samplers  lie  be¬ 
tween  the  continuous  inputs  and  D'(s),  the  inputs  can  be  separated  from 
O' (s).  The  Z  transform  of  Eq  (13)  is  given  by  Eq  (14). 

X(s)  -  Z[rf(s)]  [Xs(s)  -  Xi(»)]  ♦  Xj(s)  (14) 

Let 

D'(s)  -  Z[D'(s)] 

■  'CSS,  3 

.  I'1 

<l-s-1)(l-e-T/T  *-l) 

Since  s*1  is  a  delay  operator,  use  of  D* (s)  in  the  complementary 
filter  will  cause  the  current  estimate  of  the  actuator  position  to  be 
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based  on  its  past  position.  This  daisy  is  a  result  of  including  the 
zero  order  hold  in  0* <s)  to  provide  continuous  inputs  to  D(s).  The 
filter  will  be  implemented  in  the  digital  computer  so  it  is  possible 
to  eliminate  this  delay  by  multiplying  by  z. 

D(t)  -  z  O' (s) 


or 


(15) 


D(z)  is  not  the  Z  transform  of  D(s)  but  is  its  fioicticnal  equivalent 
for  use  In  the  digital  computer. 

Letting  K3  -  1  -  e"T^T,  the  filter  output  is  given  by 


Z(z)  - 


*3 


1  -  «3  -  1)  z-1 


(X<z)  -  Xj(z)J  ♦  Xj(r) 


The  difference  equation  for  the  filter  is  then 


X(k)  -  Z3  [Xs<k)  -  Xx(k)]  ♦  (l-K3)[x(k-i)  -  Xx(k-1)]  ♦  Xx(k) 

(16) 

It  is  interesting  to  note  that  Eq(16)  can  be  shown  to  be  a  constant 
gain  optimal  observer  for  the  gun  servo  subsystem. 
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Ill  Analysis  of  the  Gun  Servo  Subsystem 


> 


The  gun  servo  subsystem  frequencies  are  all  above  half  the  sampling 
frequency.  Sue  to  Shannon's  sampling  theorem  (Ref  8)  these  frequencies 
are  misrepresented  in  the  Z  domain;  the  effect  of  these  frequencies 
will  appear  to  be  at  frequencies  between  0  and  10  hertz  for  the  20  hertz 
sample  rate.  As  a  result,  the  gun  servo  subsystem  analysis  must  be  per¬ 
formed  in  the  continuous  time  domain. 

In  analysis  of  the  gun  servo  subsystem,  the  Delco  design  will  be 
taken  as  a  baseline  and  changes  to  this  baseline  will  be  examined.  The 
baseline  response  of  Figure  9  us as  the  given  GSS  with  a  Kql  of  100  and 
including  the  differential  pressure  compensation.  The  digital  gain  is 
20  and  the  digital  rate  feedforward  is  included.  Ho  sensor  filters  are 
included  in  the  baseline.  Notice  that  the  system  azimuth  and  elevation 
commands  are  indicated  by  as  asterisk. 

In  the  following  sections  the  effects  of  the  differential  pressure 
compensator,  feedforward  compensator,  and  rate  feedforward  will  be 
examined.  The  analyses  are  based  an  S  plane  root  locus  methods  and 
simulation  results. 

piffarentlal  Pressure  Compensator 

The  differential  pressure  (AP)  loop  is  intended  to  provide  damping  of 
the  high  frequency  structural  modes.  This  is  Important,  for  as  Figure  5 
shows,  the  83  hertz  structural  mode  is  dominant  at  the  sensor  location. 

The  AP  compensator  is  evaluated  using  the  section  of  the  gun  servo 
subsystem  shown  by  Figure  10.  Using  block  diagram  reduction,  it  can 

i 

i 
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EL  ( MRflO )  flZ  ( MRflO ) 
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*ni*5*i<  v 


SENSOR  RESPONSE 


rfoo  o!  40  o!  80  IT20  rTeo  2!  00 
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Figure  9*  Baseline  Gun  System  Response* 


Figure  10.  Differential  Pressure  Compensation  Loop. 


be  shown 


Gp(s)  m  ftia,)v  »  — 
^  R_(s)  1  ♦ 


g2  g3 _ 

►  (I-G3)  G2Hp 


Figure  11  shows  the  root  locus  of  (I-G3)  G2Hp  •  -1.  The  primary 
effect  of  the  compensator  is  to  pull  the  high  frequency  structural  poles 
to  the  left»  increasing  their  damping.  Also*  the  frequency  of  the  servo¬ 
valve  poles  is  reduced  as  they  are  moved  to  the  right. 

Xhe  Gp(s)  for  a  gain  (K*p)  of  3.15  x  1CT3  is  given  below. 


Gp(s)  - 


_ 5.17  x  1011  (s2  ♦  22.6s  ♦  5690) (s  ♦  500) (s  ♦  2500) _ 

S(s2e22.3s44940)(s2+667se296000)(s2e568se765000)(s+l500)(s422i0) 

(18) 


The  AP  compensator  has  increased  damping  of  the  high  frequency  mode  and 
slightly  reduced  the  first  mode  frequency. 
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Figure  li.  Differential  Pressure  Compensator  Root  Locus 


t 


Figures  12a  and  b  compare  the  time  response  of  the  gin  and  actuator 
with  and  without  the  compensator*  Figures  12c  and  d  make  the  same  com* 
pariscn  for  the  system  response.  Table  II  gives  rise  time  (Tr),  settling 
time  (Ts)»  peak  time  (Tp)  and  peak  overshoot  (Mp)  for  an  impulse  to 
Gp(s)  for  both  the  two  mode  and  four  mode  dynamics.  Table  III  contains 
the  same  information  for  a  step  input  cornu  and  to  the  gun  simulation. 


Table  II 

Differential  Pressure  Compensator  Gun  Response* 


Model 

Tr 

TP 

Ts 

Mp 

Without  AP 

2  Mode 

0.00206 

0.00533 

0.131 

1.50 

4  Mode 

0.00223 

0.00628 

0.308 

1.71 

With  AP 

2  Mode 

0.00246 

0.00617 

0.149 

1.23 

4  Mode 

0.00260 

0.00720 

0.405 

1.38 

*Unit  step 

input 

Table  III 

Differential  Pressure  Compensator  System  Response* 

Tr 

TP 

*s 

Mp 

With  AP 

0.05 

0.054 

0.158 

26.99 

Without  AP 

0.05 

0.054 

0.108 

26.73 

*25  arad  step  input 
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1.26- 


Figure  12.  Differential  Pressure  Loop  Time  Responses 
a.  Sensor  Response  with  no  Compensator. 


Sensor  Response  vith  Compensator 


EL  ( MRRD)  RZ  (HRRD) 

-SO. 00  -2S.00  0.00  25.00  50.00  D  -50.00  -25.00  0,00  25.00  50.00 


l!~00  oTTo  CL20  0T3Q  0T4O 

TIME  (SEC) 


EL  ( MRflD )  RZ  CMRRD) 

-50.00  -25.00  0.00  25.00  50.00  0  -50.00  -25.00  0.00  25.00  50.00 


SENSOR  RESPONSE 


1.00  0.10  0.20  0.30  0 


TIME  (SEC) 


Although  it  appears  from  liable  III  that  response  has  been  degraded 
with  the  compensator*  this  is  due  to  the  simulation  using  no  sensor 
filters  for  this  portion  of  testing.  Figures  12a*  b*  c  and  d  demon¬ 
strate  that  the  high  frequency  oscillations  are  effectively  damped. 


Feedforward  Compensator 

The  feedforward  compensator  is  an  optional  compensator  to  provide 
additional  low  frequency  gain  and  structural  damping  if  ground  firing 
tests  indicate  they  are  necessary  (Ref  3).  In  the  present  design*  how¬ 
ever,  Gj (s)  is  simply  a  gain*  KqL. 

Using  the  Gp(s)  obtained  in  the  previous  section*  one  gun  servo 
subsystem  channel  can  be  represented  by  Figure  13(a).  Figure  13(b) 
shows  the  system  after  block  diagram  manipulation.  The  closed  loop 
transfer  function  is  therefore  given  by 


Ggss<S> 


X(s)  _  (s  +  G1)GP 
^g(s)  s(l  ♦  GjGp) 


(19) 


It  is  interesting  to  note  that  if  Gp(s)  is  taken  to  be  —  (i.e.* 
neglect  structural  and  valve  dynamics)*  then  Ggfis(s)  is  given  by 


1 

s 


(  s  ♦  G1 
S  ♦  Gi 


) 


so  G^(s)  has  no  effect  on  the  rigid  body  dynamics. 

Returning  to  the  true  Gp(s)  and  letting  G^(s)  -  Kqj,*  then 


Ggss<*>' 


(s  +  Kql)  Gp 
s(l  ♦  Kql  Gp) 


(20) 
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Figure  13.  Gun  Servo  Subsystem  Reduction. 

The  root  locus  of  Kq^  Gp  ■  -1  is  given  by  Figure  14.  As  the  gain  is 
increased,  the  dominant  poles  are  no  longer  associated  with  the  high 
frequency  structural  mode  but  become  associated  with  the  servovalve. 

The  limit  on  gain  for  stability  is  Kol  ■  390.  Table  IV  gives  the  steady 
state  tracking  error  of  the  system  for  a  ramp  input.  Different  values 
of  open  loop  gain  and  rate  feed  forward  (RFF)  are  presented.  Table  V 
gives  the  simulation  time  response  figures  of  merit  for  various  open 
loop  gains. 

The  settling  time  is  very  sensitive  to  small  changes  in  the  high 
frequency  components  of  the  gun  response  since  no  sensor  filters  are 
employed.  Realizing  this.  Tables  IV  and  V  indicate  that  the  open  loop 
gain,  Kql,  has  very  little  effect  on  the  system  response. 


Table  IV 

Gun  Servo  Subsystem  Ramp  Errors 


Kol 

With  RFF 

%  Error 

No  RFF 

0 

0.23 

CO 

10 

0 

10 

100 

0 

1 

200 

0 

0.05 

Effect  of  Kql 

Table  V 

on  Gun  System  Performance* 

Kol 

Tr 

TP 

Ts 

Mp 

0 

0.052 

0.074 

0.208 

26.77 

10 

0.052 

0.072 

0.306 

26.84 

100 

0.050 

0.054 

0.158 

26.99 

200 

0.052 

0.056 

0.058 

26.56 

*25  mrad  step  input 


The  closed  loop  transfer  function  for  the  gun  servo  subsystem  with 
Kql  -  100  Is 


_ 5,17  x  1QU(s2+22.6s  ♦  5690)(s+500)(s»2500)(s+100) _ 

<s2*17 . 78s+5390)  (s2+460s«.255, 000)  (s2+166s+1, 320,000)  (s+104)  (s<*>1513  )  (s+2200) 

(21) 


33 


Effect  of  Rate  Commands  and  Rate  Feedforward 

Hie  reason  for  rate  conmand  inputs  to  the  gun  servo  subsystem  is 
shown  by  Figure  15.  This  figure  shows  the  response  of  the  gun  servo 
subsystem  without  the  integrator  (g'  (s)  -  6(i  (s))  to  a  step  input. 

gSS  o**™ 

Although  the  overshoot  would  not  be  quite  so  high  in  the  actual  system 
due  to  rate  and  acceleration  limits,  it  is  still  an  undesirable  feature. 
The  large  overshoot  of  the  gun  servo  subsystem  is  reduced  by  giving 
ramp  position  commands  consisting  of  integrated  step  rate  commands  from 
the  digital  computer. 

The  effect  of  feeding  this  rate  command  into  the  actuator  input  is 
demonstrated  by  Tables  IV  and  VI.  Table  IV  shows  that  the  steady  state 
ramp  following  error  of  the  gun  servo  subsystem  is  sero  when  rate  feed¬ 
forward  is  used.  Table  VI  indicates  that  the  rate  feedforward  improves 
the  gun  transient  response. 


Table  VI 

Effect  of  Rate  Feedforward  on  Gun  Response* 


Tr 

TP 

Ts 

Without  RFF 

0.066 

0.108 

0.168 

28.6 

With  RFF 

0.05 

0.148 

0.158 

26.7 

*25  mrad  step  input 
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IV  Analysis  of  Digital  Portion  of  Gun  System 


The  gun  servo  subsystem  discussed  in  the  last  section  receives 
rate  commands  from  the  mission  computer.  This  section  will  discuss  the 
means  of  generating  these  commands.  The  effects  of  various  elements  of 
the  digital  portion  of  the  gun  system  and  interface  elements  will  be 
evaluated  using  Z  domain  root  locus  methods  and  the  digital  simulation. 
For  the  analytical  portion  of  the  analysis,  one  channel  of  the  gun 
system  is  modeled  as  shown  in  Figure  16a. 


Figure  16.  Discrete  Gun  System  Model. 

System  Gain 

The  effect  of  the  system  gain  Kj  on  system  response  was  examined 
using  a  discrete  root  locus.  The  Z  transforms  were  found  using  the 
impulse  invariance  transformation  function  of  the  design  program  TOTAL. 
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lhe  program  uses  an  algorithm  which  first  obtains  the  partial  fraction 
•xpansion  of  F(s),  then  finds 


F(s) 


n  Al  t 

Z  — - 

to-1  *  -  e“8kT 


(22) 


where  sk  is  the  s  domain  pole  and  is  its  associated  partial  frac¬ 
tion  coefficient. 

Z  Transform  of  Gun  Servo  Subsystem.  The  gun  servo  subsystem 
including  the  zero  order  hold  is  given  by 

Gg(s)  ■  gzch*s*  Ggss^8^ 

where 

l  -  »"*T 

GzohCs)  -  - - J -  (23) 

and  T  is  the  sample  time  of  0.05  seconds. 

lhe  discrete  system  can  then  be  modeled  by  Figure  16b  where 

Gg(s)  -  Z  [GZ0H(s)  Ggss(a)] 

but  since  z  •  e*T 

Gg(s)  -  (i-s“l)  Z  [  |*GgSg(s)]  (24) 


lhe  s  transform  of  A  Ggss(s)  cannot  be  found  directly  using  TOTAL 
since  the  algorithm  of  Eq  (22)  cannot  handle  repeated  poles.  The 

sz 
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term  which  Is  formed  by  the  zero  order  hold  and  integrator  poles  is 
expanded  out  as  shown  below. 


i  «.„<•)  -  i.  -  UsiL&il 


s 


2  (1  ♦  Gp(s)  KoL 


-  ♦ 


N(s) 


®2  s(l  ♦  Gp(s)  Kqi.) 


(24) 


where 


A  ■» 


(.  ♦  Kql)  G„<«)  __  l 

1  ♦  GpCs)  Iqi 


s-0 


Substituting  A  -  l  in  Eq  (24)  and  solving  for  M(s), 


H(s)  -  (s  Gp(s)  -  l)/s 


Thus 


Gg(s)  -  Z  [  L. 


G„(s)  -  1 


s2  s2(l  ♦  Gp(s)  KoL) 


-  fcic-SL.  *«4 — )] 

t  (s-1)2  s2(l  ♦  Gp(s)  Kol) 


(25) 


The  Z  transfora  of  [s  Gp(s)  -  i]  /  (1  ♦  Gp(s)  K^)  was  fowd  using 
TOTAL.  However,  numerical  problems  were  encountered  due  to  the  extremely 
small  poles  of  Gg(z).  This  problem  was  overcame  by  a  short  program 
using  the  IMSL  subroutine  ZPOLR  (Ref  9)  to  find  the  roots  of  the  cam* 
piste  system.  The  discrete  gun  servo  subsystem  is  given  by  Eq  (26). 
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0.0509  (z2+  1.11*  +  0.416) (z2-  1.53  x  10“8  *  +  7.65  x  10"16) 

(*  ♦  0.0217X*  -  0.007 13 )(*  ♦  1.92  x  10*8) 

G  (*)  -  _ to  -  1.43  X  10-33) (g  .  lt37  x  10-48)  (26) 

(*2+  1.17*  ♦  ,445)(*2+  1.77  x  10-5*  +  l.Oi  x  10-10) 

(*2  ♦  1.25  x  10-7  z  +  4.09  x  10"15)  (*  -  .00544) 

U  -  1.43  x  10-33) (Z  .  l#37  x  10-48) (*  .  1) 

Sine*  *  ■  0  corresponds  to  i  »  -  «,  nost  of  the  terms  In  Eq  (26) 
ere  negligible  so 

•  05 (*2+  1.11*  ♦  0.416) (*  ♦  0.0217)(*  -  0.00713)(*  ♦  1.92  x  10“8) 

G  (*)  »  ■—  ■■ . .  ‘  ... . . . 

to-l)to2+  1.17*  ♦  0.445)(z2+  1.77x10-5+1. 01xl0-10)(a.0.00544) 

(27) 

The  discrete  root  locus  for  the  g*  servo  subsystem  Is  shown  In 
Figure  17.  The  maximum  gain  for  stability  (  |(z  ||  <  1)  Is  36.  For  a 
gain  of  20,  the  closed  loop  transfer  function  Is  given  by 

Xto)  _  (g«0.00713)to»0.0217)(*2e  1.11*  ♦  0.416) (*  ♦  1.92  x  10*8) _ 

(*2+0.0458z  ♦  0.0183) (*2+  1.24*  ♦  0.460)(*+  1.54xl0”8) (*-0.00777) 

(28) 

The  response  of  this  model  Is  compared  against  the  results  of  the 
simulation  In  Table  VII. 

Table  VIII  gives  the  dominant  pole  locations  for  several  values  of 
gain.  Notice  that  the  first  mode  poles  are  still  essentially  cancelled 
by  the  zeros  at  -0.55  ±  j  0.33  so  have  little  effect.  The  time  response 
corresponding  to  these  system  poles  can  be  determined  by  referring  to 
Appendix  E. 


39 


Figure  17.  Discrete  Gun  System  Root  Locus 


Table  VII 

Comparison  of  Analytical  and  Simulation  Response* 

<KX  -  20) 


k 

6<k) 

Simulation 

(mrad) 

Analytical 

0 

0 

0 

1 

25.4 

25.4 

2 

23.8 

23.8 

3 

26.5 

26.2 

4 

23.9 

24.2 

5 

25.4 

25.4 

*25  mrad  step  input 

Table  VIII 

Poles  of  Closed  Loop  Gun  System  as  a  Function  of  Gain 

*1 

Dominant  Poles 

10 

0.49 

-0.60  ±  j 

0.31 

12.5 

0.36 

-0.60  ±  j 

0.31 

15 

0.21 

-0.61  db  j 

0.30 

17.5 

0.080  ±  j 

0.096 

-0.61  ±  j 

0.29 

20 

0.022  ±  j 

0.13 

*-0»  62  db  j 

0.28 

22.5 

-0.036  db  j 

0.14 

—0.62  ±  j 

0.26 

Additional  insight  into  the  effect  of  Kj  is  obtained  by  completely 

neglecting  the  gun  servo  subsystem  dynamics.  The  Z  transform  of  the 

gin  servo  subsystem  is  then  simply  .  The  root  locus  is  then  a 

z- 1 

line  to  the  left  from  z  -  1  crossing  z  -  0  at  Kj»*  20  and  z  »  -1  at 

40.  This  would  Indicate  instability  at  a  gain  of  40  and  a  deadbeat 
response  for  a  gain  of  20.  The  simulation  indicates  instability  at  a 
gain  of  39,  so  in  this  respect  the  simple  model  is  better  than  the  more 
complex  model  of  £q  (27).  Ihis  is  probably  due  to  numerical  errors  in 
finding  Gg(z)  and  deletion  of  the  small  roots. 


Digital  Rate  Feedforward 

This  section  will  discuss  the  need  for  digital  rate  feedforward  in 
the  mission  computer.  Figure  18  shows  the  response  of  the  system  with 
a  gain  of  20  to  a  ramp  input  without  rate  feedforward.  The  gun  output 
lags  the  command  input  by  the  sample  time  (0.05  seconds). 

There  are  two  methods  of  dealing  with  this  problem.  Hie  first  is 
to  give  the  gun  position  commands  which  are  one  sample  period  in  advance; 
however,  this  places  an  additional  burden  on  the  target  prediction 
algorithms.  The  other  alternative  is  to  use  the  digital  rate  feedfor- 
ward  as  discussed  in  section  II. 

Use  of  the  rate  feedforward  x'educes  ramp  following  errors  from  5 X 
to  zero. 

Sensor  Filters 

Up  to  this  point,  none  of  the  simulations  have  used  a  filter  on  the 
sensor.  Although  the  system  works  well  without  a  filter,  the  system 
performance  can  be  Improved  through  their  use. 

Analog  Sensor  Filter.  As  a  result  of  Shannon's  sampling  theorem 
(Ref  8),  any  sampled  input  having  a  frequency  greater  than  one  half  the 
sampling  frequency  (fs)  will  appear  to  have  a  different  frequency.  In 
fact,  it  will  have  a  frequency  between  0  and  fs/2.  This  effect  is  called 
aliasing  and  in  the  case  of  the  gun  system,  increases  the  settling  time 
by  causing  low  frequency  (  <  fs/2  )  oscillations. 

To  prevent  this,  a  low  pass  filter  can  be  placed  on  the  input  to 
the  mission  computer.  Reference  3  suggests  a  first  order  lag  filter 
with  a  10  Hz  cutoff.  When  such  a  filter  is  used,  the  system  overshoot 
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Figure  18.  Response  of  Gun  System  with  No  Digital  Rate  Feedforward. 
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and  settling  time  Is  increased  as  shown  in  Figure  19.  This  degradation 
o£  performance  is  due  to  decreasing  the  bandwidth  of  the  system. 

The  output  of  the  sensor  filter  with  a  cutoff  of  10  Hz  is  down 
3.7  dB  at  11.6  Hz,  thus  the  first  mode  is  not  attenuated  much  more  than 
the  desired  10  Hz  signal  which  is  down  3  dB.  Fortunately,  the  sensor 
does  not  sense  much  of  the  first  or  second  modes  so  the  sensor  filter 
cutoff  frequency  can  be  increased.  As  shown  in  Table  IK,  a  sensor 
filter  of  30  Hz  provides  better  performance  than  the  10  Hz  filter. 


Table  IX 

Effect  of  Sensor  Filters  on  System  Response* 


Cutoff 

Frequency 

Tr 

TP 

Ts 

Mp 

None 

0.05 

0.054 

0.158 

26.49 

10  Hz 

0.05 

0.104 

0.268 

32.32 

30  Hz 

0.05 

0.136 

0.156 

27.06 

*25  mrad  step  input 


Hie  information  in  Table  IX  indicates  that  either  sensor  filter 
cutoff  frequency  Increases  the  peak  overshoot.  This  is  because  the 
servo  bandwidth  is  reduced  from  100  Hz  to  the  sensor  cutoff  frequency. 

A  smoothing  filter  on  the  output  of  the  zero  order  hold  will  have  the 
same  effect,  thus  its  cutoff  frequency  should  also  be  somewhat  greater 
than  the  10  Hz  suggested  by  Reference  3. 

Digital  Sensor  Filter.  The  outputs  of  the  analog  sensor  filter 
are  Inputs  to  the  digital  filtering  algorithm.  This  algorithm  estimates 
the  position  of  the  gun  based  on  the  measurement  and  the  integrator  out¬ 
put  which  provides  rigid  body  response. 
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EL  ( HRflO )  RZ  ( MRRD ) 

-SO. OO  -25.00  0.00  25.00  50.00  0  -SO. 00  -25.00  0.00  25.00  50.00 


Figure  19.  Effect  of  10  Has  Sensor  Filter  on  System  Response. 
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The  difference  equation  for  the  digital  filter  (Eq  (16))  was 
employed  in  the  simulation  program*  Table  X  compares  the  system  perfor¬ 
mance  for  different  filter  gains.  For  a  gain  (K3)  of  0  the  filter  com¬ 
pletely  ignores  measured  outputs.  With  a  gain  of  1,  the  measured  outputs 
are  used  directly. 


Table  X 

Effect  of  Digital  Filter  Gains  on  System  Response* 
(4%  settling  criteria) 


*3 

Tr 

Ts 

“p 

0 

0.05 

0.054 

0.058 

26.99 

0.5 

0.05 

0.054 

0.100 

26.99 

1 

0.05 

0.054 

0.454 

26.99 

*25  mrad  step  input 


A  4%  settling  criterion  was  used  since  a  5%  criterion  indicated  no 
difference  in  the  system  performance  for  different  gains.  This  would 
indicate  that  the  errors  caused  by  aliasing  are  less  than  5%.  Although 
the  best  response  is  obtained  for  K3  -  0,  this  gain  could  not  be  used  in 
the  actual  system  since  the  system  would  be  running  open  loop.  The 
best  gain  for  use  in  the  system  must  be  selected  based  an  noise,  and 
system  error  considerations. 

Computation  Time 

When  implemented  on  the  mission  computer,  output  measurements  and 
control  Inputs  corresponding  to  a  given  sample  period  are  separated  by 
some  finite  time  delay.  This  computation  time  (Tc)  is  used  for  analog 
to  digital  and  digital  to  analog  conversions  as  well  as  time  to  perform 
digital  computations. 
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When  this  time  delay  is  included  in  the  system,  the  single  channel 
model  is  given  by  Figure  20a,  where  C  is  the  computation  time.  A  delay 
in  the  computer  output  is  the  same  as  a  delay  in  the  system  output  so 
the  system  can  also  be  given  by  Figure  20b. 


Figure  20.  Gun  System  Model  Including  Computational  .^alay. 


If  Ggfs(s)  is  allowed  to  be  l/s  (Ignoring  complex  dynamics),  then 
Gg(s)  -  G2QH  (  i  )  e  C® 


(  )(  j  )(e“Cs) 


(29) 


1-e  _-Cs 

x  e 


Using  the  modified  Z  transform  (Ref  1:234)  the  delay  can  be  treated 
by  letting 
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Z  [x(t-C)]  -  *m*  Z  [X(s)  a^Ts  ] 


X(z,m) 


where  a  is  one  plus  the  integer  number  of  sampling  periods  delayed  and 
A  is  a  number  between  0  and  1  such  that 


C  -  (m  -  A)  T 


For  our  case  it  is  reasonable  to  consider  computation  delays  no  greater 
than  one  time  period  so  m  ■  1  and 

A  -  1  *  *  <30 

T 


Z[x(t-c)]  -  Zmlz  l  e^Ts  ]  x^fc) 


-  t*1  <l-s’1)  Z  [  ]  Xcd(z) 

-1  -1.  AT  ♦  T  (l-A)  zml  _  .  . 

-  *  <1-* ” a,a-i)2 - 


-  at  [z  »  (1/A  -  1)' 

»<*•!) 


*cd<1> 


The  complete  forward  transfer  function  is  then  given  by 

X(z)  m  Ki  XA  [»  ♦  (1/A  -  D] 

COO  *(z-l) 


where  A  i«  given  by  Eq  (31). 


For  Kj  m  20  and  T  ■  0. 05f 


the  closed  loop  response  is  given  by 


> 


X  (z)  Az  ♦  (1  -  A) 

y— 7  t  ■  . . . .  (34) 

XCd(z)  z2  „  z  +  (1-A) 

Table  XI  compares  the  output  of  the  simulation  with  the  prediction 
of  Eq  (34).  The  theoretical  results  match  the  simulation  results  except 
for  the  transients  caused  by  neglected  dynamics. 

Figure  20  shows  the  gun  response  for  a  computation  time  of  0.024 
seconds.  Table  XII  gives  the  figures  of  merit  for  the  simulation  with 
different  computation  times.  Both  show  that  computation  time  produces 
increased  overshoot  and  settling  time.  Results  indicate  that  the  system 
gain  must  be  reduced  to  compensate  for  computation  times  in  excess  of 
0.005  seconds  (10%  of  the  sampling  time). 


Table  XI 

Comparison  of  Simulation  and  Analytical  Computation  Time  Effects* 


k 

Tc 

9  00 

•  .01 

Tc 

-  .024 

Analytical 

Simulation 

Analytical 

Simulation 

0 

0 

0 

0 

0 

1 

20,0 

18.1 

13.0 

9.4 

2 

29.0 

28.2 

31.2 

30.2 

3 

26.8 

28.0 

33.8 

35.0 

4 

29.6 

23.9 

26.2 

26.7 

5 

24,6 

25.0 

21.4 

21.5 

*25  mrad  step  input 
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EL  CMRRO)  RZ  ( MRRD) 

-SO. OO  -25.00  0.00  25.00  50.00  _  -SO. 00  -25.00  0.00  25.00  50.00 


Figure  21.  Effect  of  Computation  Time  on  System  Response  (T.-.024  sec). 


50 


Table  XII 

System  Response  for  Different  Computation  Times* 

<Kl  -  20) 


c 

Tr 

TP 

TS 

Mp 

0 

0.05 

0.054 

0.158 

26.99 

0.006 

0.056 

0.144 

0.162 

26.91 

0.01 

0.064 

0.114 

0.218 

29.67 

0.024 

0.082 

0. 128 

0.432 

37.75 

*25  mrad  step  input 


Word  Length  Considerations 

Due  to  a  finite  word  length,  a  digital  control  system  must  repre¬ 
sent  its  data  in  quantitized  units.  The  effects  of  vordlength  in  the 
computer  are  highly  dependent  on  the  specific  algorithms  used,  so  will 
not  be  discussed  further.  However,  the  word  length  of  the  data  inter¬ 
face  can  be  checked  to  determine  its  accuracy. 

Although  the  data  busses  can  carry  16  bits  of  information  per 
channel,  the  A/D  and  D/A  converters  have  a  10  bit  word  length  (Ref  3), 
The  ratio  of  the  largest  number  to  the  smallest  which  can  be  repre¬ 
sented  by  the  word  is  given  by 


!aa 2S.  • 

valn 


2n 


where  n  is  the  number  of  magnitude  bits.  Subtracting  one  sign  bit  from 
the  10  bit  word,  9  bits  are  available  for  magnitude  information. 


m  29  m  512 

^min 

If  the  largest  angle  to  be  represented  is  52  mrad,  then  the  minima 
resolution  is  52/512  or  .1  mrad.  nils  is  well  within  the  requirement  of 
the  ,5  mrad  static  accuracy  specification. 
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V 


Analysis  of  Muzzle  Response 

Up  to  this  point,  the  gun  system  response  has  been  examined  only 
at  the  sensor  location.  This  section  will  deal  with  the  response  of 
the  gun  at  the  muzzle.  The  muzzle  response  will  be  examined  and  a 
controller  designed  to  improve  damping  of  the  first  structural  mode. 

Effect  of  Structural  Modes  on  Muzzle  Response 

By  examining  Figure  22  which  shows  the  gun's  structural  mode 
shapes,  it  can  be  seen  that  the  first  (11.6  Kz)  mode  has  a  very  small 
displacement  at  the  rear  of  the  gun  and  a  low  slope  at  the  pivot  posi- 
tion.  This  indicates  that  neither  the  actuator  position  sensor  nor  the 
angle  resolver  at  the  pivot  can  sense  the  first  mode  response.  This  is 
fortunate  as  simulation  results  indicated  unacceptable  oscillations  for 
larger  first  mode  output. 

The  gun  muzzle,  however,  lies  at  the  point  of  maximum  mode  shape 
slope.  It  is  the  muzzle  which  determines  the  projectile  direction  so 
the  first  mode  oscillation  has  a  significant  effect  an  gun  performance. 
The  muzzle  frequency  response  of  Figure  23  also  shows  the  high  gain  of 
the  first  structural  mode  at  the  muzzle. 

Figure  24  shows  the  gun  muzzle  time  response  which  can  be  compared 
against  the  sensor  response  of  Figure  9.  Two  results  can  be  seen.  For 
a  step  input  the  muzzle  has  large  overshoot  and  icng  settling  time;  for 
a  sinusoidal  input  the  muzzle  tracks  the  input  quite  well.  This  means 
that  the  first  structural  mode  will  cause  significant  oscillation  in 
target  acquisition  but  target  tracking  will  not  be  seriously  affected. 
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It  may  be  acceptable  to  have  a  two  second  settling  period  during 
target  acquisition.  If  not,  seme  form  of  compensation  is  required  to 
increase  the  damping  of  the  first  mode.  The  next  section  will  describe 
the  preliminary  design  of  such  a  compensator. 

Controller  for  Structural  Modes 

This  preliminary  structural  mode  compensator  design  makes  two 
important  assumptions.  First  is  that  the  controller  is  necessary  and 
second  is  that  the  barrel  cluster  is  non-rotating.  The  rotating  barrel 
cluster  will  cause  cross  coupling  between  the  axes  due  to  gyroscopic 
precession.  Although  it  is  not  known  what  the  magnitude  of  this  effect 
will  be,  it  is  suspected  that  it  may  be  significant. 

Three  different  controller  coif igurat ions  were  considered,  a  digital 
controller,  use  of  the  existing  feedforward  compensation  and  a  feedback 
compensator. 

The  structural  modes  cannot  be  controlled  by  the  digital  controller 
for,  as  mentioned  earlier,  the  natural  frequencies  of  the  gun  lie  outside 
the  controllable  region  of  10  Hz.  An  alternative  approach  would  be  to 
consider  an  intermediate  digital  processor  which  would  operate  at  a 
higher  frequency.  A  controller,  having  a  sample  rate  of  100  Hz,  was 
designed  to  control  the  first  structural  mode,  but  aliasing  of  the  third 
and  fourth  inodes  caused  unacceptable  response. 

As  shown  in  Chapter  III,  the  feedforward  compensator,  Gj(s),  has 
little  effect  cn  the  low  frequency  response  of  the  gun  so  this  compen¬ 
sator  will  not  be  capable  of  controlling  the  first  mode.  Thus,  feedback 
compensation  becomes  the  only  viable  method  of  first  mode  control. 
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The  sensors  cannot  directly  measure  the  muzzle  position  since  they 
have  very  low  amplitude  response  to  the  first  mode.  It  is  conceivable 
that  a  high  gain  bandpass  filter  could  pick  out  the  first  mode  response 
from  the  sensor?  but  such  an  approach  would  require  an  extremely  high 
order  filter.  Ideally?  a  sensor  for  first  mode  control  should  be  lo¬ 
cated  at  the  muzzle.  The  rotating  barrels?  however?  make  this  a  very 
complex  hardware  problem?  not  to  mention  the  difficulties  of  working 
with  a  rotating  coordinate  frame.  It  is  apparent  then  that  some  form  of 
observer  is  required  to  find  the  barrel  response  based  on  sensor  output. 

A  conventional  observer  for  the  Y  channel  of  the  gun  servo  sub¬ 
system  has  the  form  (Ref  6) 

i  -  [A]  £  ♦  B  Yc  ♦  L  [Ys  -  CS  X  ] 

(35) 

ym  -  cm£ 

where  A  Is  the  system  matrix,  E  is  the  input  matrix,  L  is  the  observer 
gain  matrix  and  and  C§  are  the  servo  and  muzzle  output  matrices. 

A  A 

X  is  the  estimated  state  vector  and  Ys  and  Y^  are  the  sensor  output  and 
estimated  muzzle  output?  respectively.  This  was  considered?  but  it 
appears  that  at  least  a  five  or  possibly  seven  state  system  model  would 
be  required  to  include  the  integrator  and  gun  structural  modes. 

Hie  method  chosen  to  determine  muzzle  response  was  to  realize  that 

&<!>  -  i  ,  V - -  (36) 

Xc(s)  s  (s2  +  2f  Wls  +(i)j  ) 
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neglecting  higher  frequency  dynamics,  end  that  the  sensor  transfer 
function 


*s<s>  *  I 

XcCs) 

Eqs  (36)  and  (37)  yield 

Xs(s)  s2  +  2 1  ♦  (t)^2 


(37) 


(38) 


Feedback  of  the  estimated  muzzle  position  causes  the  system  roots 
to  move  toward  the  right  half -plane  as  the  gain  is  increased,  so  the 
rate-acceleration  feedback  of  Eq  (39)  was  employed. 


Hs(s) 


s(g  ♦  10) 
s2  ♦  2{  (OjS  +  djj2 


(39) 


It  would  be  desirable  to  place  this  type  of  compensator  in  an  inner 
feedback  loop,  perhaps  the  same  as  the  AP  compensator,  leaving  a  unity 
feedback  outer  loop.  However,  for  convenience,  the  compensator  was 
placed  in  the  outer  loop  as  shown  by  Figure  25. 

The  root  locus  of  Figure  26  shows  that  as  open  loop  gain  is  in¬ 
creased,  the  dominant  first  mode  poles  move  to  the  left.  A  of  0.003 
gives  a  |  of  0,17.  Increasing  the  gain  further  increased  the  damping 
ratio,  but  overall  system  performance  was  degraded  due  to  reducing  the 
system  bandwidth.. 
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Using  KoL  m  0.003,  gun  servo  subsystem  transfer  function  is 


Ggss<s> 


0.454  (s2  ♦  7.29s  ■»  5312' 
s<s2  ♦  24.7s  ♦  5066) 


after  elimination  of  non-dominant  poles  and  zeros. 

Finding  the  Z  transform  of  Eq  (40)  ir.  series  with  a  zero  order 
hold,  as  was  done  in  section  IV, 


Gg(z)  -  0.04_28  (z  ♦  1.25)  f  0.426  (41 

(z-1)  (z2  +  Iz  ♦  0.290) 

Die  root  locus  of  Figure  27  indicates  that  the  system  damping 
ratio  for  a  system  gain  of  20  is  about  0.25  which  is  somewhat  low,  A 
gain  of  15  produced  the  time  response  shown  in  Figure  28.  This  figure 
and  Table  XIII  show  that  the  compensator  does  help  the  step  response 
of  the  muzzle  by  reducing  settling  time  and  peak  overshoot.  However, 
the  square  root  of  the  mean  square  muzzle  error  (MSE)  for  a  2  Hz  sine 
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Figure  26.  Structural  Mode  Compensator  Root  Locus 


EL  ( MRflO )  FIZ  (MRfiD) 

-50.00  -25.00  0.00  25.00  50.00  n  -50.00  -25.00  0.00  25.00  50.00 


Figure  28a.  Structural  Mode  Compensator  Time  Response. 
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SENSOR  RESPONSE 


sc* 


input  is  increased  by  the  compensator,  indicating  that  the  tracking 
ability  of  the  gun  has  been  reduced. 

The  structural  inodes  of  the  gun  can  be  controlled  as  demonstrated 
by  this  low  level  compensator  design. 


Table  XIII 

Structural  Mode  Compensator  Response* 


Ts 

Tr 

Tp 

Mp 

t]  MSE  ** 

Servo  Responses 

Without  Compensator 

0.158 

0.05 

0.054 

26.99 

- 

With  Compensator 

0.09 

0.102 

0.104 

25.91 

- 

Muzzle  Response 

Without  Compensator 

>2 

0.048 

0.07 

39.28 

3.62 

With  Compensator 

0.156 

0.070 

0.150 

26.45 

4.69 

*25  mrad  step  input  **20  mrad  2Hz  sine  input 

Additional  simulations  indicate  that  10g  errors  in  the  estimation  of 
first  mode  frequency  and  damping  ratios  have  little  effect  on  the  con¬ 
troller’s  ability  to  damp  step  inputs.  The  primary  drawbacks  of  this 
preliminary  filter  design  are  that  low  frequency  tracking  capabilities 
have  been  degraded  and  barrel  cluster  rotation  has  been  ignored. 

It  is  possible  that  a  closed  position  loop  could  improve  this 
response.  In  addition,  a  lag-lead  compensator  could  be  added  at  Gj  to 
allow  a  higher  KQL.  Barrel  rotation,  however,  would  be  a  much  more 
complex  problem  to  deal  with. 

The  difficulties  associated  with  implementing  a  realistic  compen¬ 
sator  for  the  first  structural  mode  appear  to  be  more  than  the  problem 
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warrants .  Since  the  tracking  response  Is  acceptable  with  no  compensa¬ 
tion,  It  nay  be  best  to  simply  wait  for  the  acquisition  transients  to 
die  out  before  firing  the  gun. 
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VI 


The  major  results  of  this  study  will  now  be  summarized  and  recom¬ 
mendations  made  based  on  these  results. 

Summary 

The  gun  system  easily  meets  all  specifications  developed  in  the 
Introduction.  The  gun  system  settles  in  under  0.2  seconds  and  has  an 
overshoot  of  less  than  10%  for  a  system  gain  of  20,  independent  of  Kq^ 
and  with  or  without  a  differential  pressure  compensator.  Although  the 
full  dynamics  of  the  gun  servo  subsystem  are  complex,  it  can  be  modeled 
as  a  pure  integrator  for  discrete  analysis.  The  slew  rate  and  acceler¬ 
ation  limits  have  little  or  no  effect  on  the  linear  analysis  at  the 
gains  used. 

A  system  gain  of  20  provides  essentially  a  deadbeat  response  for 
the  gun  system.  As  computation  time  exceeds  10%  of  the  sampling  time 
it  is  necessary  to  reduce  the  system  gain  to  stay  within  overshoot 
specif ications.  In  order  to  meet  ramp  following  specifications,  a  rate 
feedforward  must  be  Included  in  the  mission  computer. 

Ihe  digital  sensor  filter  has  little  effect  on  the  gun  response 
using  the  specifications  in  the  Introduction;  however,  if  response  is 
examined  at  a  closer  level,  the  filter  does  reduce  the  system  settling 
time.  Use  of  a  10  Hz  analog  filter  for  the  sensor  output  severely 
degrades  performance;  however,  a  30  Hz  filter  eliminates  high  frequency 
Inputs  to  the  discrete  controller  while  maintaining  a  reasonable  time 
response. 
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The  first  gun  structural  mode  causes  overshoots  of  over  50%  in  the 
muzzle  angle  and  settling  times  in  excess  of  2  seconds  for  a  step  input. 
The  effect  of  the  first  mode  is  not  as  severe  in  tracking  polynomial  or 
sinusoidal  inputs.  If  barrel  cluster  rotation  is  neglected,  a  compen¬ 
sator  can  be  designed  to  control  the  acquisition  (step  input)  oscilla¬ 
tions.  This  controller  degrades  tracking  of  other  inputs,  however. 

The  complexity  of  dealing  with  barrel  cluster  rotation  and  the  uncom¬ 
pensated  tracking  response  indicate  that  the  acquisition  problem  is  best 
handled  by  allowing  the  barrel  vibrations  to  damp  out  before  the  gun  is 
fired. 

The  M61  movable  gun  is  an  extremely  fast  system  which  does  not 
require  any  compensation  provided  that  commands  to  the  gun  servo  sub¬ 
system  are  rate  commands  as  in  the  current  configuration.  The  addition 
of  differential  pressure  compensation,  sensor  filters  and  selection  of 
a  good  KqL  all  improve  the  gun  system  response  by  small  degrees,  result¬ 
ing  in  a  highly  effective  gun  pointing  system. 

Recommendations 

This  study  has  shown  the  M61  movable  gun  control  system  to  be 
effective,  neglecting  barrel  cluster  rotation.  It  appears  that  addi¬ 
tional  investigation  into  the  combined  effects  of  first  mode  oscillation 
and  barrel  cluster  rotation  is  required.  Investigations  in  this  area 
could  also  lead  to  a  realistic  method  of  reducing  acquisition  settling 
time. 
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Appendix  A 


Derivation  oL  Gun  Dynamics  Model 


In  the  following  discussion,  the  gun  dynamics  model  will  be  derived 
based  on  a  lumped  parameter  model  which  might  be  used  in  a  finite  ele¬ 
ment  analysis. 

The  equation  of  motion  for  a  lumped  parameter  model  of  a  viscously 
damped  structure  is  (Ref  16:390) 

[M]  *£  ♦  [C]  g  +  [X]  q  -  Q  (42) 


where  M,  C,  and  K  are  the  mass,  damping  and  stiffness  matrices,  respec¬ 
tively,  q  is  the  generalized  coordinate  vector  and  Q  the  generalized 
force  vector. 

The  mass  matrix  is  symmetric  and  positive  definite  so  it  can  be 
expressed 


M  .  MhL»t 


Letting  q 


M‘%  i , 


'£  ♦  M  i  *  M  s  - 


M'4  9 


where 

M  -  W'H  M  W'% 
M  -  W4  M  [mT* 


(43) 
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It  can  be  shown  that  if  A  can  be  expressed 

M  -  2  2  ttrp  [B]P/  (44) 

r-1  p-0 

where  is  a  different  coefficient  for  each  combination  of  r  and  p, 

then  the  system  can  be  decoupled  into  its  normal  modes.  The  transform 

maticn  to  decouple  the  system  is  given  by  u  -  v  where  f  is  the 

transformation  matrix  composed  of  the  eigenvectors  of  the  system  and  has 
•1  T 

the  property  •  ■  •  , 

Making  the  transformation,  the  system  is  given  by 

£  +  [2*  w]  u  ♦  [d)2]a  -  [»]T  [Mr*  Q  (45) 

where  (2^  (si]and[(i)2Jare  both  diagonal  matrices  If  Eq  (44)  is  satisfied. 

One  mode  of  Eq  (45)  is  given  by 

ut  4  »  (dj  in  *(t)i2  uj  -  tj1  [m]"1  a  -  X^(Pa)  f  (46) 


where  is  the  eigenvector  associated  with  the  1th  mode  and  Xi(Pe)  is 
defined  as  the  mode  shape  at  the  actuator  where  the  force  f  is  applied. 
Equation  (44)  is  satisfied  if  l  ^  is  constant  for  i  -  l»n. 

The  system  output,  y,  can  be  given  by 


y(t) 


n 

2  X.  (Ps)  ut(t) 
i-i 


(47) 


where  Xi(Ps)  is  the  mode  shape  evaluated  at  the  output  location.  Taking 
the  Laplace  transform  of  Eq  (46)  and  employing  Eq  (47), 
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Determination  of  Gun  Output  Coefficients 


The  gun  structural  dynamics  are  based  on.  the  frequency  response 
plots  provided  in  Reference  5*  Two  sets  of  dynamics  are  derived;  one 
is  for  the  muzzle  location,  and  the  other  for  the  sensor  location. 


Mizzle  Dynamic  Equation 

From  Appendix  A,  the  gun  muzzle  output  equation  can  be  given  by 


%(a) 


*M<S> 

TOT- 


n 

2 

i-t 


Alm 

(s2  ♦  2{  (j)iS  +  Qjj2) 


(49) 


where  Aj^  are  the  muzzle  output  coefficients  to  be  determined,  {  is  the 
gun  damping  ratio  to  be  determined  and  Wj  are  the  first  four  normal  fre- 
quencies  given  by  Reference  5  as 


(i)£  ■  11.6  Hz 

(i)2  -  74.9  Hz 
(1)3  ■  68.7  Hz 
<1)4  -  83.5  Hz 


72.88  rad/sec 
219.2  rad/sec 

431.6  rad/sec 

524.6  rad/sec 


Figure  29  shows  the  frequency  response  for  the  gun  muzzle.  The 
dashed  line  indicates  the  response  given  in  Reference  5  and  the  solid 
line  the  response  of  the  model  developed  here.  The  figure  indicates 
the  maximum  response  is  located  near  11.6  Hz  and  has  a  value  of  20  dB 
corresponding  to  a  magnitude  tyg)  of  10. 
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Figure  29a«  Modeled  Muzzle  Frequency  Response 


FREQUENCY  l HERTZ  ) 


Using  the  relation  (Ref  2:299) 


% 


then 


2  „  1  "  V  1  ■  (1/Mm)2 
;  "  2 


(50) 


so 


i  -  0. 05 


We  also  know  that 


lim  G>j^(s)  -  lim  2 


s  ->  0 


Ai 


a-»  0  i-1  s2  +  (jjjs  +  ^ 


-  1 


(51) 


to  satisfy  the  steady  state  (rigid  body)  requirements.  Thus 


2  -  1  (52) 

i-1  (i)i2 

If  the  expression  for  the  magnitude  of  G3^(j  (i))  were  written  as  a 
function  of  {  ,  (j)i»  (i)»  and  Aj,  it  would  be  possible  to  take  three 
measurements  of  (j)  and  Mq  from  the  frequency  response  and  form  three 
equations  in  four  unknowns  (A^ , i  -  1,4).  If  £q  (52)  were  added,  four 
equations  in  four  unknowns  could  be  solved  for  the  A^'s.  This  method 
was  tried  for  finding  A}  but  the  nonlinear  algebraic  equatims  became 
unwieldy,  so  a  more  conventional  synthesis  method  was  used. 

The  gun  transfer  function  is  given  in  its  factored  form  as 

3  2  2  \ 

TT  (sZ  ♦  2tj  (d^S 

0^(8)  -  K  ^ -  (53) 

IT  (s2  ♦  2(  (sUS  ♦  &h2) 


where 


the  same  as  the  gun  damping  ratios,  (j)^  ,  (i^  »  and  (i)^  were  then 
chosen  such  that  the  frequency  response  of  Gm(s)  was  fit  to  the  true 
gun  response,  as  shown  by  the  solid  line  of  Figure  29.  Note  that  the 
magnitude  is  not  well  matched  at  high  frequencies.  In  order  to  match 
the  magnitude  response,  zeros  would  be  required  near  20  and  45  Hz; 
however,  the  phase  diagram  does  not  show  zeros  at  these  locations. 

The  partial  fraction  expansion  of  Eq  (53)  was  then  found  using  the 
design  program  TOTAL.  The  form  of  this  expansion  is 


4 

TT 

i-1 


Aj  +  B[S 

s2  +  2 {  GOjS  ♦  (i)i2 


Since  B  does  not  affect  the  steady  state  response  and  was  small 
compared  to  A,  it  was  set  to  zero.  A  check  of  the  frequency  response 
Indicated  that  no  change  from  the  factored  form  was  detectable* 

The  muzzle  output  coefficients  are  then  given  by 
Aj  ■  6041 
A2  —  -6605 
A3  -  -1319 
A*  -  1918 
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Sensor  Dynamic  Equations 

Two  different  sensors  were  considered  by  Delco.  The  first,  a 
linear  sensor  is  mounted  on  the  actuator.  The  second,  an  angle  sensor, 
is  located  an  the  gun  gimbals.  The  frequency  response  of  the  two 
sensors  is  nearly  identical,  so  one  transfer  function  is  considered  for 
both  sensors. 

Hie  sensor  output  coefficients  were  determined  using  the  same  pro- 
cedure  as  for  the  muzzle.  Figure  30  shows  the  given  and  modeled  fre- 
quency  responses. 

The  sensor  output  coefficients  are  given  by 
A}  «  356.2 
A2  -  0 

A3  -  48170 
A^  -  185624 
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Appendix  C 


Gun  Servo  Subsystem  State  Variable  Representation 

the  gun  servo  subsystem  dynamics  must  be  put  into  state  variable 
form  for  the  computer  simulation.  A  physical  variable  form  composed  of 
phase  variable  blocks  was  chosen  to  represent  the  system. 

the  state  variable  model  is  formed  by  first  considering  each  trans¬ 
fer  function  individually  and  then  connecting  inputs  and  outputs  accord¬ 
ing  to  the  functional  block  diagram  of  Figure  31. 


Figure  31.  State  Variable  Block  Diagram. 
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a  Variable  Form  (Ref  13) 

Hie  phase  variable  form  places  the  transfer  function 
U(s)  bft  s11  ♦  .  .  .  b2  s2  +  bt  s  +  b0 


into  the  form 


Y(s)  sJ*  ♦  .  .  ,  a2  s2  ♦  ai  s  +  ac 


xx  -  x2 


x2  -  X3 


*n  “  “aoxl  "  alx2  •  •  •  *  u(t) 

y(t)  -  (b0  -  +  (bL  -  a1bn)x2 

*  *  •  *  <bn-l  “  ^i-l  bn)xn  +  bn  u(t) 


'««  ■  ^ 


Using  Eq  (54) 


X1  "  U1 

yi  -  xi 


G2(s)  - 


Y2(s)  _ Ay 

S3  ♦  2/  (|)yS^  ♦  (t)y2S  ♦  0 


Using  Eq  (54) 


0  0  1  X3  +  0  u2 

0  -d)y2  -2{{j>y  X4  1 

J  L  J  L  . 


n  -  Dv 


0] 


*2 

x3 

*4 


(56) 


Differential  Pressure  Compensator 


Hp(s)  - 


Y3(s) 


U3(s)  (s  ♦  1/TjXs  +  1/t2) 


Since  the  differential  pressure  compensator  has  no  complex  eigen¬ 
values,  it  is  easily  implemented  in  a  Jordan  form  (Ref  18). 

Taking  the  partial  fraction  expansion  of  Eq  (56), 


H  (s) 
P 


AcC 


A  ,  B  “I 
1  ♦  1/Ti  1  +  1/T2  ' 


where 


A  -  T2/(T2  -  TX) 
B  -  Tj/OTj  -  T2) 


The  state  variable  representation  is  then 


x5 

x6 


1/Ti  0 
o  i/T2 


x5 

Ac 

♦ 

J 

> 

Ac 

u3 


x5 

x6 


y3  -  [  A  B  ] 

Dote  that  the  simulation  program  uses  T  *  1/r 


(57) 


Gun,.  Dynamics 

The  Gun  Dynamics  are  treated  as  a  set  of  decoupled  modes  each 
represented  by  a  phase  variable  block. 
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G3(s)  -  «  £ 


U4CS)  l-l  <s2  ♦  2}  (ijjs  +  o)^ ) 
Che  of  these  phase  variable  blocks  is  given  by 


*8 


0  1 
-dJj2 


" 

x7 

♦ 

0 

x8 

1 

-  M 

(58a) 


Similar  blocks  exist  for  states  9  through  16.  The  output  equation  is 
given  by 


y u  “  A1  x7  +  h.  x9  +  ^  X11  +  A4  *13  +  A5  x15  <58b* 

where  the  output  coefficients,  A^,  will  be  different  for  sensor  and 
muzzle  output. 


Feedback  Compensator 


H  (s)  -  Ys(s)  -  Cfls^  C  f2g2  *  cf3s  ♦  Cf4 

8  U5(s)  s2  +  Cf5s2  ♦  Cf6s  +  Cf7 


Using  Eq  (54) 


'  • 

' 

x17 

• 

x18 

■ 

• 

x19j 

_ 

0  10 
0  0  1 
-Cf7  -Cf6  -Cf5 


“ 

'  * 

x17 

0 

X18 

♦ 

0 

J 

x19 

.  _ 

1 

:  j 

y5  "  Ccf4*  cf7Cf 1  cf3“  Cf 6cf 1  cf2”  cf5cfll 


x17 

x18 

x19 


(59) 


♦  Cflu5 
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■oral 


Feedforward  C 


Bagasse 


Gi(s) 


Zft  . 

»6 


JJL 


♦  C12  s  ♦ 


c13 


s  *  Cj4  s  ♦  C15 


^20 

X21 


u6 


y6 


Cc12  -  C15  Cu 


c12  *  c14  cll] 


*20 

X21 


Sensor  Filter 


Gf(s)  IZ.  C2 1  s2  ♦  c22  s  ♦  C23 

°7  s2  ♦  C24  s  ♦  C25 


*22 

*23 


0  1 
-c25  -C24 


*22 

0 

♦ 

*23 

1 

"7 


y7  "  Cc23  *  C25  C21  C22  “  c24  c2 1^ 

From  Figure  31  it  can  be  seen  that 
Uj  ■  u 

«2  ■  u  *  y6  *  y3 

u3  "  y4s  "  y2 


x22 

x23 


u4  "  y2 

u5  "  y4s 

«6  "  *01.  ^1  *  y5^ 


♦  C11  u6 


*  C2l  u7 


4 


and 


^sensor  "  y4 
ymuzr  le  "  y4a 
?f  liter  “  y7 


where  y4  uses  the  sensor  output  coefficients  in  Sq  <5 8b)  and  y^  uses 
the  muzzle  output  coefficients. 

Using  Eqs  (55),  (56),  (57),  (58),  (59),  (60),  (61),  and  (62),  the 
entire  state  variable  representation  is  given  by 


*2 

*3 

*4 


*  lvrl/T  *6  +  °11  K°L  Cf  1 
2  1 


(Aj  x7  4  A2  x9  4  A3  Xjj 


♦  44  xX7  ♦  A5  xl5) 


♦  C11  *0L  C(Cf4  “  cf5  Cfl>  x17  *  <cf3  *  Cf6  CfP  x18 


♦  (Cf2  -  Cf5  Cfl)  xl9] 

4  (C13  -  Cl5  Cu)  X2(J  4  (CI2  -  CI4  Cjj)  X21 

*5  "  ml^l  x5+  Ac  C“V2  ♦  a1x7  ♦  a2x9  *  ^ll  ♦  A4X13  ♦  Vis! 
X6  -  -l/Tl  X6  4  Ac  t-\X2+  AlX74  A2X94  A3Xu4  A^*  AjXjj] 


86 


» 


*8  " 


-  2J(i)lxg 


x9  -  X10 

X10  "  Avx2  “  w22  x9  "  ^^2  x10 
X11  "  X12 

XI2  -  V2  "  ^32  X11  "  2^3  x12 
x13  "  x14 

*14  "  V2  "  V  x13  "  2H  x14 
x15  "  X16 

x16  “  V2  "  w52  x15  “  2*u5  x16 
x17  "  x18 

x18  "  x19 

x19  **  ^ix7  *  ^3X9  ♦  AjX^  ♦  A^Xj3  ♦ 

*  Cf7  X17  "  Cf6  x18  *  Cf5  x19 
x20  "  *21 

*21  "  *0L  X1  "  Cfl  ^1X7  *  ^2X9  +  ^11  *  A4x13  +  ^*15* 

-*0i[<cf4-cf7cfl)x17  *  (Cf3’Cf6CCl)x18  *  (Cf2"Cf5Cf  1*19^ 
-c15x2o  -  C14x21 
x22  -  *23 

*23  "  Alx7  *  *2X9  *  A3x11  +  A4X13  *  A5X15  "  C25x22  *  C24x23 

The  output  equations  are  then  given  by 

ysensor  "  CAlx7  *  *2X9  ♦  ^11  ♦  Vl3  ♦  VlJ 
y barrel  -  l>blx7  ♦  ^2*9  ♦  ^,3*11  *  *b4*l3  *  ^5*15^ 
yfilter  "  <c23  *  *22  ♦  <C22  *  htPll*  *23 

♦  CjiCAjX?  ♦  AjXj  ♦  AjXjj  ♦  A^Xjj  ♦  A^x^) 
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Appendix  D 
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M61  Simulation  Program 

The  M61  Simulation  Program  is  designed  to  provide  gun  system 
response  as  a  function  of  time  for  a  variety  of  commanded  Inputs.  The 
program  is  coded  in  CDC  FORTRAN  IV  and  requires  64K  core  memory  for 
execution.  The  simulation  requires  approximately  IS  seconds  of  CP  time 
per  second  of  simulation  time. 

The  main  program  is  responsible  for  data  input ,  output  and  reduc¬ 
tion,  and  acts  as  the  simulation  executive.  All  gun,  control,  and 
simulation  parameters  are  read  from  the  input  file  in  free  format.  The 
main  program  then  echos  the  inputs  and  processes  them  for  use  in  the 
simulation.  After  initializing  values,  the  simulation  begins. 

The  program  uses  the  integration  package  ODE  (Ref  8)  to  integrate 
the  GSS  state  equation  between  sampling  instants.  The  actuator  outputs 
are  found  at  the  sensor,  barrel,  and  sensor  filter  locations  and  con¬ 
verted  to  gun  angles  by  subroutine  LNTQAN.  At  the  sample  instant,  the 
control  input  is  formed  by  subroutine  CONTROL.  If  T,.  is  not  zero,  the 
input  is  delayed  by  the  appropriate  number  of  time  intervals.  Note  that 
the  computation  time  may  not  be  exactly  as  specified  due  to  the  quanti¬ 
sation  of  OT,  the  integration  interval.  DT  is  selected  to  provide  2S0 
data  points  for  the  CALCGMP  plotter  based  on  the  final  time. 

After  the  simulation  is  complete,  the  azimuth  and  elevation  errors 
for  both  the  sensor  and  barrel  location  are  computed  and  printed  along 
with  the  commands  and  positions. 
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Rise  time,  settling  time,  peak  time  and  peak  overshoot  are  deter* 
mined  for  the  azimuth.  Mean  and  mean  squared  errors  are  calculated  for 
elevation.  Subroutine  PLOTTER  then  plots  the  barrel  and  sensor  angles 
if  desired. 

Subroutine  GUN  contains  the  state  equation  for  the  X  and  Y  channels 
of  the  gun.  Since  the  channels  are  identical,  the  same  equations  are 
used  for  each  channel.  Actuator  position,  rate,  and  acceleration  limits 
are  incorporated  in  the  servovalve  equations. 

Subroutine  LNTQAN  transforms  the  actuator  length  to  gun  angles  for 
obtaining  GSS  outputs.  Subroutine  ANTOXN  transforms  gun  angles  to 
actuator  lengths  for  the  digital  controller.  This  subroutine  returns 
both  actual  length  and  change  in  length  from  center. 

Subroutine  COMMAND  generates  position  and  rate  commands  for  the  gun 
system.  The  azimuth  and  elevation  can  be  commanded  independently  with 
several  different  functions,  including  steps,  ramps,  parabolas,  sines, 
cosines,  and  a  random  acceleration  input.  Subroutine  NOISE  generates 
the  random  input. 

Subroutine  PLOTTER  plots  gun  azimuth  and  elevation  positions  and 
position  comnands  against  time  using  standard  CALCCMP  routines. 

Subroutine  CONTROL  uses  the  algorithm  of  section  II  in  simulating 
the  mission  computer.  Kg  is  used  as  a  switch  for  the  rate  feedforward. 
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PROGRAM  M61A1  (INPUT, OUTPUT. PLOT.  RESULT, TAPE5=RESULT) 
OIMENSI3M  T (253) ,AZM(253) ,ELM <253 ) , A7C < 253) ,ELC<253> 
DIMENSION  EL  (  253)  ,  A7  ( 253) 
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PRINT*,-  “ 

PRINT*, "CONTROL  VARIABLES  ARE* 
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XLI=X(  1) 

YLI=X( 2k) 

30  CONTINUE 
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Z  and  S  Plane  Relationships 

Pole  locations  in  the  Z  plane  can  be  related  to  S  plane  poles  and 
associated  response  using  Figure  32.  Lines  of  constant  a  in  the 
plane  correspond  to  circles  centered  at  the  origin  in  the  Z  plane  and 
lines  of  constant  damped  frequency  are  radical  lines.  While  stability 
in  the  S  plane  is  indicated  by  left  half  plane  poles,  Z  domain  stability 
is  indicated  by  poles  with  a  magnitude  less  than  1. 


Figure  32.  Z  Domain  Unit  Circle 
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